function iirTest()
  % Load the reference data
  data = load("gse2.txt");
  data = double(data);
  % Set an IIR filter
  b = [0.000401587491686, ...
       0.0, ...
      -0.001606349966746, ...
       0.0, ...
       0.002409524950119, ...
       0.0, ...
      -0.001606349966746, ...
       0.0, ...
       0.000401587491686]';
  a = [1.000000000000000, ...
      -7.185226122700763, ...
       22.615376628798678, ...
      -40.733465892344896, ...
       45.926605646620146, ...
      -33.196326377161412, ...
       15.023103545324197, ...
      -3.891997997268024, ...
       0.441930568732716]';
  b2 = [0.005028154401050966,
        0.0,
       -0.010056308802101932,
        0.0,
        0.005028154401050966];
  a2 = [1.0,
       -3.787291923503399,
        5.384566634300907,
       -3.407019620837378,
        0.8097462844356603];
  % Filter data
  H = dfilt.df2(b, a);
  H2 = dfilt.df2(b2, a2);
  yfilt = H.filter(data);
  yfilt2 = H2.filter(data);
  %yfilt = filter(b, a, data);
  % Write the data
  fileID = fopen('iirReference1.txt','w');
  for i=1:length(yfilt)
     fprintf(fileID, '%.13e\n', yfilt(i));
  end
  fclose(fileID);
  fileID = fopen('iirReference2.txt','w');
  for i=1:length(yfilt)
     fprintf(fileID, '%.13e\n', yfilt2(i));
  end
  fclose(fileID);
end
